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ABSTRACT 

One of the most fundamental features of a black hole in general relativity is 
its event horizon: a boundary from which nothing can escape. There has been 
a recent surge of interest in the nature of these event horizons and their local 
neighbourhoods. In an attempt to resolve black hole information paradox(es), 
and more generally, to better understand the path towards quantum gravity, 
“hrewalls” have been proposed as an alternative to black hole event horizons. In 
this paper, we explore the phenomenological implications of black holes possessing 
a surface or “hrewall”, and predict a potentially detectable signature of these 
hrewalls in the form of a high energy astrophysical neutrino hux. We compute the 
spectrum of this neutrino hux in diherent models and show that it is a possible 
candidate for the source of the PeV neutrinos recently detected by IceCube. 

This opens up a new area of research, bridging the non-perturbative physics of 
quantum gravity with the observational black hole and high energy astrophysics. 

Subject headings: Phenomenology of quantum gravity. Physics of black holes. 
Quantum aspects of black holes. Radiative transfer; scattering. Neutrino, muon, 
pion, and other elementary particles; cosmic rays. Accretion and accretion disks 


1. Introduction 

One of the most celebrated predictions of General Relativity has been the possibility of 
forming black holes: spacetime singularities that are surrounded by event horizons. These 
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event horizons are boundaries in the spacetime from which nothing can escape. In contrast 
to black holes, other astrophysical compact objects such as neutron stars, possess a physical 
surface that shows visible signs of radiation to the outside universe. In this work we explore 
the possibility that black holes, too, may not be so invisible after all. 


There has been much recent interest in studying the nature of black hole event horizons 
and their local neighbourhoods, in an attempt to resolve black hole information paradox(es), 
and more generally to better understand the path towards quantum gravity (e.g., Almheiri 


et ah 2013 Vaz| 2015 Braunstein, Pirandola fc Zyczkowski 2013; Hutchinson fc Stojkovic 


2013t [Pen &: Broderick]|2014t |Page|[20T^ [Lowe &: Thorlacius||2013t [Banks fc Fischlei1|2013 


Lee & Yeom 2013 Kawai, Matsuo &: Yokokura 

2013 Giddings 2013[ [Hsu 2013 Kim, Lee 

& Yeom 2013 Harlow & Hayden 2013 Susskind 

2013; Jacobson [2014 ; [Larjo, Lowe & Thor- 


lacius 2013 ; [Giddings 2013 Nomura fc Varela|2013 Rama 2014} Hwang, Lee fc Yeom 2013 


Nomura, Varela fc Weinberg 2013; Avery, Chowdhury fc Puhm 20131 Hossenfelder 2012 


Susskind]|2012t |Brustein]|2012t [Giveon fc Itzhaki||2012t [Bena, Puhm fc Vercnocke||2012t jOri 


2012t [Banks fc Fischler]|2012t [Susskind][2012t [Chowdhury fc Puhm[[2012[ [Mathur fc Turton 


2014 Nomura, Varela &Weinberg 2013 Bousso 2012 Saravani, Afshordi & Mann 2015). In 


particular, alternatives to classical black hole event horizons such as firewalls (Almheiri et 
al. 2013), fuzzballs, (e.g., Mathur & Turton 2014), or other exotic byproducts of quantum 


gravitational collapse (e.g., Giddings][2014 Dodelson fc Silverstein[[2015 ) have been discussed 
in the literature. It is generally believed that, in lieu of an event horizon, accretion onto 
(or scattering off) such exotic byproducts may have astrophysical observable consequences 
that distinguish them from classical black holes ([Broderick, Loeb fc Narayan 2009 Pen & 


Broderick 

2014, 

Broderick et al. 

2015 


We know that a quantum theory of gravity is non-local on the scale of the Planck 
length/time, if not longer, which means that a radical change in the structure of space-time 
within a Planck length of the horizon can affect the microscopic vicinity outside the event 
horizon. In this work, we assume that gravitational collapse leads to a “hrewall” (a term 
we shall use as a placeholder for any exotic alternative to classical event horizons), that can 
show visible signs of emission to the outside universe. Such emission may be possible because 
an accreting firewall is not in global thermal equilibrium, as the timescale for global thermal 
equilibrium is generically much longer than the timescale for local radiative processes for 
macroscopic objects. This is analogous to meteorites hitting the Earth’s atmosphere, where 
light is emitted from the meteor (within seconds), while the long timescale of heating up the 
entire earth and global thermal equilibrium are irrelevant. 


We construct a phenomenological model to describe and constrain the nature of this 
radiation and its spectrum. A key feature of our prediction for this radiation is that it will 
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primarily consist of neutrinos. This is because, as in core-collapse supernovae, which nom¬ 
inally share similar radiative properties to hrewalls (e.g. surface density and composition), 
neutrino propagation has a much smaller optical depth and therefore transports energy out 
more efficiently. 


An interesting implication is that accretion onto black holes/firewalls could be a possible 


source for the recently detected high energy (PeV) neutrinos by IceCube (Aartsen et ah 


2014). The source of this detection is as yet unknown, and several candidate sources have 
already been ruled out or strongly constrained. For example, hadronuclear pp processes 
(the mechanism of neutrino production in intergalactic shocks, starburst galaxies, and some 
active galactic nuclei (AGN) models) require power laws harder [p < 2.1 — 2.2) than those 
favored by IceCube’s detection (Murase, Ahlers & Lacki 2013). Inner jets of blazars may be 
able to account for IceCube’s PeV signal, but they have difficulty explaining the sub-PeV 
neutrino events (Murase et ah 2014). Further constraints on the candidate sources are set 
by arguments that no extragalactic object that we know of could function as the origin of 
both the PeV neutrinos and the ultra high energy cosmic rays (Yoshida & Takami 2014[). 


We start by outlining the main features and assumptions in our model in Sec. We 
then summarize our estimate of the total flux emitted by the black holes (Sec. 1^ followed by 
modelling their spectra (Sec. |^. While we consider both blackbody and power law spectra, 
we subsequently focus on the power law spectrum in Sec. (which are more interesting at 
higher energies), comparing our results to observations and other models. Finally, Sec. 
concludes the paper. 


2. Modelling Accretion onto a firewall 


The black holes in our model are accretion powered, have Keplerian flows, and have 
reached approximate steady state. 

We assume that the accreting black holes must show evidence of surface radiation. 


Broderick, Loeb & Narayan (2009) considered the possibility of photon emission from the 
surface of Sagittarius A*. As their predicted photon flux exceeds the Sag. A* flux in mil¬ 
limetre and infrared observations, they conclude that there must be an event horizon. In 
this work, we argue that at the energy scales involved, the neutrino cross-section is much 
smaller than the photon cross-section, hence (similar to core-collapse supernovae) most of 
the energy dissipated in an optically thin region will be carried out by the neutrinos. Based 
on this assumption, we go on to estimate the total flux of neutrinos from all the black holes 
(or hrewalls) in the universe. 
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We estimate that a fraction eg (where 0 < < 1) of the gravitational binding energy 

of particles falling onto the snrface of the black hole is radiated away (and reaches inhnity), 
while the rest goes into growing the mass of the black hole. For a continnons accretion flow, 
as we are assnming here, this gives a total Inminosity as measnred at infinity of: 


Loo = egMc^. 


( 1 ) 


A fraction of the total luminosity, L^, is converted by the accretion disc into electro¬ 
magnetic radiation: 

^Loo- (2) 

The remainder of Loo, will be eventually radiated at the surface as neutrinos. We call this 
surface luminosity as measured at infinity L^-. 

Ly = Loo ~ = e^Mc^, (3) 


where we have introduced a new efficiency e^ = eg — egf] for the surface radiation that we 
will be interested in for this work. In terms of the observed total integrated flux: 

Mc2 




(4) 


where di is the luminosity distance. Now we want to sum over (|^ to find the total flux 
contribution from all the black holes in the observable universe. We treat the contributions 
from the stellar mass black holes and supermassive black holes separately. We estimate 
that 3.3% (Fukugita & Peebles 2004) of the stellar mass budget goes into forming stellar 
mass black holes, while Mbh ~ 10 Mq and Msmbh ~ 10® Mq are typical masses of stellar 
mass and supermassive black holes respectively. As a typical accretion rate we take M ~ 
0.1 MEdd = 0.1 X 16LEddC“^, where LEdd is the Eddington limit for the black hole mass. 


3. Total Flux 


Let us make the following definitions of characteristic radii 

GM 


rh = 


(^1 V'l - , 


Oi — 1 -|- (1 — a^)^'^®((l -|- -f (1 — 

«2 = ^3a2 al, 

GM r . 

Hsco — —^ 3 -|- 02 ~ v (0 ~ tti)(3 -|- Oi -|- 202 ) 


(5) 

( 6 ) 

(7) 

( 8 ) 
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where is the location of the event horizon, and risco is the radius of the innermost stable 
circular orbit. For simplicity, we assume that the black holes in our model are typically 
spinning with a dimensionless spin parameter of a* ~ 0.7. 


For the disc efficiency we use Frank, King & Raine (2002): 


1 / 2M 

V Srisco 


(9) 


For example, for eg = 0.5, we will have efficiencies of r; = 0.2 and e^ = 0.4 for both 
stellar mass and supermassive black holes. 


We primarily use models from Springel & Hernquist (2003) and Di Matteo et ah (2003), 


for the evolution of the star formation rate (SFR) and the supermassive black hole accretion 
rate (BHAR) comoving densities, although our final results are not sensitive to our choice 
of these models. Both SFR and BHAR are fitted by a double exponential function: 


p{z) = PO] 


ga(2 Zm) 


+ a ’ 


( 10 ) 


where the best fit parameters for the SFR are a = 3/5, b = 14/15, Zm = 5.4, po = 
0.15 MQyr~^Mpc~^, and for the BHAR they are a = 5/4, b = 3/2, = 4.8, po = 

3 X 10“"^ MQyr~^Mpc~^. 

The total flux then is 


e,.c 


^ / (0.033 psFRc") 


O Jo 


X 


47r(l + + (1 + z^^VLm 


dz, 


( 11 ) 


pSMBHr 1 _ 
^ ujot L'^maxJ — 


X 


eyC 




o Jo 


Pbharc 

1 


47r(l + zy^/ff^+~(l^^^)WiM 


dz, 


( 12 ) 


where we use Ho = 73kms ^Mpc Ha = 0.7, and VLm = 0.3. We are working in comoving 
coordinates, so the volume elements in (pT| and (|I^ are 


dV 


c (1 + zfdl 


dz, 


(13) 
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where d-A = „ ^ ^ ^ dz and di = (1 + z^dA- 

Ho{l+z) JO ^QA+(l+^)3nAf ^ ^ 

For Zjriax = 9, near the redshift of reionization, Eq’s (p!I| and (12) yield nentrino flnxes 


< 9] 

F^:^r[z < 9 ] 




rs> 


3.8 X 10 erg s ^cm ^sr 

8.9 X 10“® erg s“^cni“^sr“^. 


(14) 

(15) 


4. Spectrum 

To compare to observations, we need information abont the spectrnm of the nentrino 
flnx. We consider two different models: a power law spectrnm and a blackbody spectrnm. 


4.1. Power Law with p ~ 2 — 3 

For the power law model, the observed nentrino phase space density ftot[Ei,] is 
, , r~ (p(r.)c‘dV\ { irr\ \ (E,{1 + z)' 

= J (xM^rj 

where the valne of is fixed by the Inminosity: 


(16) 


L„(p) = 4irr^ 


dEp j^3 


-p-2 


^ 1 dEp, 


(17) 


TT^C^h^ \Eo 

ta is the apparent radius of the black hole, and Ep = 1.96 x 10^® erg is the Planck energy 
(even though the result is insensitive to this choice for p > 2). Note that, given that 
the neutrino phase space density cannot exceed unity due to the Pauli exclusion principle, 
E Eo is the natural lower cut-off for a power-law spectrum. Details of how to obtain va 
can be found in Broderick, Loeb Sz Narayan (2009): for a* ~ 0.7 we have ta ~ 5GMc~^. 


The results depend on the precise choice of index p, which is one of our parameters. We 
will focus on p ~ 2 — 3. The number flux is 


dN 


dt dA dEy dVt 

where 5 is a normalization constant set by 

Ey dN 
JEo 


= Bftot[Ey 


El 




(18) 


dt dA dE, da 1 ^““'- 


(19) 
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We make the simplifying assumption that the neutrinos emitted from the set of stellar 
mass and supermassive black holes each have approximately the same minimum energy 
and thus the same spectrum. 


An example of a power law with p ~ 2 is hrst order Fermi acceleration (Fermi 1949). 
First order Fermi acceleration is a mechanism that describes the acceleration of charged 
particles, such as electrons, crossing strong shocks. Upon meeting the shock, due to the 
disturbance in the magnetic held caused by the shock, there is a probability for the incident 
electrons to get bumped back the way they travelled. There is also a probability for them to 
pass through the shock. For those charged particles that bounce back, this process occurs 
repeatedly until they escape. Each bounce causes the particle to gain energy and the statistics 
of this process yields a spectrum of: 


dEp Eq \ Ep J 


( 20 ) 


where p ~ 2 in the nonrelativistic case (in the relativistic case, the spectrum is still a power 
law but the index can vary depending on the details of the scattering; see e.g., Vietri 2008). 
No and Eo are the initial number and energies of the charged particles. 


Fermi acceleration is ubiquitous for charged particles in astrophysical environments, 
due to their tight coupling to the magnetic held. Even though neutrinos are not charged 
and do not interact with an electromagnetic held, they do interact with a gravitational 
held which varies signihcantly on short time/length scales close to the hrewall. Therefore, 
depending on the radiative properties of the hrewall, it is possible for neutrinos to bounce 
back and forth between the accretion how and the black hole hrewall (through gravitational 
scattering), which could yield a power law whose precise index would depend on the details of 
the propagation model. We shall next discuss a possible mechanism for gravitational Fermi 
acceleration. 


4.2. Gravitational Fermi Acceleration 

Fermi acceleration requires presence of converging magnetic mirrors., which can gener- 
ically occur in a magnetohydrodynamic (MHD) turbulent medium (2nd order), or around 
strong shock fronts (1st order). For neutral particles such as neutrinos, a similar ehect can 
happen due to gravitational interactions. In particular, relativistic neutrinos can only be 
signihcantly dehected by the gravitational held in the vicinity of black holes. However, null 
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geodesics around ordinary Kerr black holes, generically either fall into the horizon or escape 
to inhnity, and thus don’t have the opportunity to up-scatter to higher energies. 


The picture can be different if we replace the black hole horizon by a firewall. In 
particular, one possibility that gives a local description of Bekenstein-Hawking entropy for 
hrewalls excises the spacetime inside the stretched horizon, imposing boundary conditions 
(Saravani, Afshordi & Mann 2015[ ). This implies that the firewall may act as a mirror for 
test particles. As motivated in Sec. [T} we shall assume that this picture is indeed valid for 
neutrinos due to their weak interaction with matter. However, photons can be absorbed by 
the hrewall, due to their electromagnetic interactions. 


We now have our two gravitational mirrors, one provided by the gravity of the black 
hole, and the other by a(n exotic) firewall. We shall next study the evolution of neutrino 
energy due to the fluctuations of spacetime geometry, caused by e.g. matter accreting onto 
a black hole. The energy is given by: 


E = ep„ 


( 21 ) 


where = Sq would be the time-like Killing vector for a stationary metric, while is the 
four-momentum. For a time-dependent metric, the rate of energy change is given by: 


dE 




1 /^ 
2 V dagi'^dxP ) 


( 22 ) 


i.e. the relative change in particle energy is roughly proportional to the relative change in 
the metric: 

where = dx'^/dX, is defined using an affine parameter A. 


The change in particle energy can be divided up into a systematic drift, and a random 
walk. The systematic change can be calculated using adiabatic invariance, i.e. that action 
variables for bound orbits of integrable systems (away from the resonances): 


J^ = 



(24) 


remain invariant under slow changes of the Hamiltonian. Geodesics in Kerr spacetime are 
indeed integrable, as three integrals of motion (energy, azimuthal angular momentum, and 
Carter constant; Carter 1968) exist. For null geodesics, the action variables simply scale as 
J^ExR(xEx Mbh, so adiabatic invariance yields: 


(AlnF^) ~ — AIuMbh, 


( 25 ) 
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i.e. bound relativistic particles systematically redshift with the size of the black hole horizon. 

The computation of the random contribution to gravitational redshift is more involved, 
as it depends on the details of metric fluctuations induced by the turbulent accretion disc. 
However, this could simply be measured by studying geodesics of test particles in existing 
general relativistic magnetohydrodynamic (GRMHD) simulations (with live geometry) of 
accretion discs in their infalling regions. Let’s try to estimate this: 


Given that an accretion disc of radius R is turbulent on the scale of its scale-height, H, 
we estimate the relative change in energy of a particle in one orbit to be: 


A In i? ~ A In \g^ 


fll/ \ 


I H ^ ^disc 

R Mbh 


(26) 


Over several orbits, \n.E would undergo a random walk (plus the systematic drift in (25)), 
with the variance: 


((AhiE)") ~ 


fH 




^disc \ ^ 


(27) 


\R J V TLbh / Grbit 

On the other hand, the accretion rate (for a radiatively efficient accretion disc) is roughly 
given by: 

Mbh ~ ass'^disc 

where ass ~ 0.01 — 0.1 is the celebrated Shakura-Sunyaev viscosity parameter (e.g., see 

(29) 


Balbus & Hawley 1998). This yields: 


((AlnE)O 


27ra 


^disC A 1 T,r 

X — X A In Mbh- 


ss 


Ml 


BH 


Equipped with Eqs. (25) and (29), we can write a Fokker-Planck equation for the evolu¬ 
tion of the energy distribution of neutrinos, trapped between a hrewall and the gravitational 
potential barrier of an accreting black hole: 


where 


dn 


d 


-An H —B 
2 


dn 

d\nE^ 




n 

A 

B 


E,. 


dN 


dN 


dEy dhiER 


V 

Mb 


dt -Mbh 

^{{A\nEf) = 
dt 2 ttoiss 


"idisc Mbh 

Mbh Mbh 


(30) 


(31) 

(32) 


(33) 
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The equilibrium solution to the Fokker-Planck equation (30), assuming zero flux through 
(high or low energy) boundaries, gives: 


cl Inn 2 A 

'0 = 1 -= 1 -~ 1 + 

^ dlnE B 


or equivalently 


p - 1 ~ 1.3 


ass \ f rUdisc/M] 


0.01/ V 0.1 


Airass 

^disc/AfBu’ 

(34) 



mu \ 

(35) 


The requirement of a significant disc to black hole mass fraction (> 10%), to fit the 
IceCube measurement of p ~ 2.5, suggests a stellar origin, e.g., the fallback discs of gamma 
ray bursts. Alternatively, the accretion rate (28) could be significantly smaller for radiatively 
inefficient flows (e.g.. Yuan, Quataert fc Narayan||2003), but the details will be more model- 


dependent in this regime. 


Note that this calculation provides the steady-state spectrum of trapped neutrinos. 
However, we expect the spectrum of neutrinos that escape to infinity to be the same, as 
the escape probability of relativistic particles is independent of energy, and only depends on 
geometric factors. 


Alternatively, protons could be accelerated via first order Fermi mechanism, and then 
produce neutrinos through hadronic processes (e.g. Waxman & Loeb (2011)). We will con¬ 
sider indices p rsj 2 — 3. 


In general, we think that a power law with p ~ 2 is a natural choice. This is because 
high energy processes (such as possible Planck-scale physics within firewalls) favour hard 
spectra, but the requirement of convergence sets the lower limit of p = 2 (equal energy per 
decade). 


4.3. Blackbody 

If we assume that the black holes are in approximate steady state and are close to 
thermodynamic equilibrium, we can consider them as blackbodies. Again, we make the 
simplifying assumptions that the set of stellar mass and supermassive black holes each have 
approximately the same surface temperature T and thus spectrum. The temperature can be 
determined by the neutrino Stefan-Boltzmann law: 


L 


u,BB 


4:7iar\ 


rsj 


(36) 
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• Frejus v^, Daum et al. 

• Frejus Vg, Daum et al. 

• AMANDA II Abbasi et al. 

IceCube 68%, Aarsten et al. 

- IceCube 90%, Aarsten et al. 

AGN, Mannheim et al. 

GRB, Hummer et al. 

— SMBH, €^,=.4, BB, This work 

— BH, €y=A, BB, This work 

— BH, €^=.045, p=2.51, This work 
- SMBH, €^=.146, p=2.36, This work 

Relic V 


Fig. 1.— Models and observations of nentrino flnxes. The AGN and gamma ray bnrst 


(GRB) predictions are from Mannheim, Protheroe & Rachen (2001) and Hummer, Baerwald 


& Winter (2012) respectively. The IceGube detection hts are from Aartsen et al. (2014), and 


the atmospheric neutrino detections are from Frejus (Daum et al. 1995) and AMANDA-II 
( Abbasi et al.| 2010). The predictions for the neutrino fluxes from surfaces of black holes 
with a blackbody and power law spectrum are from this work. 



• AMANDA II Vfj, Abbasi et al. 

- IceCube 68%, Aarsten et al. 

- IceCube 90%, Aarsten et al. 

- AGN, Mannheim et al. 

- GRB, Hummer et al. 

- BH, €,.=.045, p=2.51. This work 

- SMBH, ey=.146, p=2.36. This work 


Fig. 2.— A close-up of Fig. ([^ within the vicinity of the IceGube detection. 
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With ey = 0.4, we find Tbh = 1-2 x 10^ K and Tsmbh 
distribution is: 

dN _ I El 

dEiy ^E^IkT _|_ 


6.7 X 10^ K. The energy 

(37) 


The number fiux is again obtained using (|18|), with the phase space density 

1 


ftoAEu] = 


p{z)c^dV\ f 7ir\ 


Lu,bb J \‘^Ed\j \ ^ 


(38) 


5. Results 


The IceCube collaboration has recently observed a fiux of PeV (10^^ eV) neutrinos 


(Aartsen et ah 2014), and the source of these high energy neutrinos is as yet unknown. 
There are only upper bounds beyond 2 PeV, and these are consistent with the power laws 
we consider in this work. 


Figure [T] shows our results for the neutrino spectra for the stellar mass black holes and 
the supermassive black holes, for the two cases of power law (with p ~ 2) and blackbody 
spectra, as discussed in the previous section (Note that. Figure shows only a sample pair 
of values that fit the data. For the complete set of parameters that fit the IceCube 
data, see Figure |^. Also included in the figure are neutrino fiux predictions from models 


of AGN ( 

Mannheim, Protheroe & Rachen 

2001 

) and CRB’s ( 

Hummer, Baerwald & Winter 

2012 

), along with a fit to IceCube’s detection from ( 

Aartsen et ah 

2014 

), and atmospheric 


neutrino detections. The atmospheric neutrino observations are from Frejus (Daum et ah 


1995) and AMANDA-II ( Abbasi et al.||2010 ). The blackbody models do not have interesting 


contributions to the high energy spectra, but it would be interesting in future work to 
compare these results to upper limits on other sources of low energy neutrinos. In this work, 
we are interested in high energy neutrinos, so we focus on the power law model results. In 
Figure we have zoomed in on our results and the other predictions, in the vicinity of the 
IceCube detection. We see the firewall neutrino signal can certainly fit the IceCube detection 
and is a possible candidate for the source of these high energy neutrinos. 


Figure]^ shows the parameter combinations for neutrino emission efficiency and spec¬ 
tral index p for which our power law model yields a good fit to the IceCube data. For 
the power law models, fits to the IceCube’s data are within the physically acceptable range 
(0 < Cg < 1) at 96% and 47% probability for stellar mass black holes and supermassive 
black holes respectively (or 92% and 37% for 0 < < 0.5). These probabilities can change 

by ±5% across different models for star formation and supermassive black hole growth. 
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P 

Fig. 3.— IceCube likelihood contours for firewall neutrino emission efficiency and spectral 
index p, for supermassive (Upper) and stellar mass (Lower) black holes. The shaded regions 
from dark to light correspond to the la, 2a, and 3a contours. These contours are obtained 
from a likelihood-ratio test using the test statistic —2AlnL = —2(lnL — In Lmax)- 
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At the la confidence level the fit parameters for stellar mass (snpermassive) black holes 
~ p = 2.26 — 2.68 (2.27 — 2.45), while at the 2a level they are 

> 10-5 (10-3-5) and p = 2.1 - 2.7 (2.1 - 2.47). 

6. Conclusions 

To summarize, we have argued that accretion onto astrophysical firewalls is likely to pre¬ 
dominantly lead to a flux of cosmological neutrinos. If this spectrum has a power-law shape, 
e.g. due to an analog of Fermi acceleration for neutrinos, then it could well fit the observed 
IceCube spectrum of high energy neutrinos at reasonable efficiencies from astrophysical black 
holes. Future work should focus on more detailed spectral modelling of neutrino emission 
from putative firewall atmospheres, as well as other possible observational smoking guns for 
this scenario. It is certainly exciting to entertain the possibility that neutrino astrophysics 
could open a new window into the physics of quantum gravity and black holes. 
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